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ABSTRACT 



Submersibles are limited in depth by the strength of the pressure 
hull. The pressure hull^ in turn, is limited in strength by the properties 
of the material used in fabrication, construction techniques, and by the 
weight of material used in construction. Once the ballast wei^t has been 
taken up in hull veight, a limit is reached in depth. To go deeper with the 
same^ payload, the hull must be enlarged or size may be held constant, the 
hull strengthened and flotation material attached to support the extra 
weight. Either of these methods, enlarging the hull or adding flotation, 
increases size greatly as depth is increased. 

This thesis is a study of total vehicle size versus depth for 
wei^t-limited deep submersibles with and without buoyancy material 
attached. 
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I. DISCUSSION 



A» Yolme-Llmlted Subnarines 

1. Conrentional submarines have operated routinely 
in the upper layers of the ocean for the better part 
of a century. These submarines, designed primarily 
as military vehicles, are large, providing the neces- 
sary internal volume required by men and equipment. 
Wel^t of hull and equipment is less than the sea 
vater displaced, thus it is necessary that these 
vessels carry ballast to attain neutral bvioyancy. 

2. To design a submarine for maximum depth, it is 
necessary that the hull vel^t be increased. This 
can be accomplished with a volxime-limlted submarine 
with no increase in size by trading ballast weight 
for hull wei^t. The depth to which a volme-limlted 
submarine can be designed with no increase in size 

is then, a function of the ballast it carried. Beyond 
that depth the weight of hull and equipment exceeds 
the wei^t of sea water displaced, and the vehicle 
becomes weight -limited requiring additional buoyant 
support to maintain neutral buoyancy. 
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B. Deep SiibmersiBles 



1. Oceanographers require small laboratories to, carry 
them to the depths of tte sea. The designer's problem 
in building submarines to satisfy the needs of ocean- 
ography is to build vehicles capable of carrying a 
scientifically usefxil payload to a scientifically 
interesting depth, visually the maximum depth possible. 
Submarines of a given size_, limited in depth by 
hull vei^t and hull materials_, have several avenues of 
approach by vhich their depth may be increased. First, 
vehicles may be made larger vith payload fixed. 

Improved pressure hull materials may be used or flota- 
tion material may be attached to the vehicle to provide 
the lift necessary to restore neutral buoyancy. The 
first and last of these methods result in large craft vhich 
are less maneuverable, more coBp>licated, and more expensive 
to build. This thesis examines and compares the penalties 
paid in achieving maximum depth by the methods Just des- 
cribed utilizing various pressure hull materials. 
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II. PROCEDURE 



A. General Approach. 

1. A number of assumptions are made concerning the 
family of vehicles to he considered^ the structiiral 
materials available, the desired payload, the pro- 
pulsion characteristics and the other features vhich 
vill affect the huoyancy-vei^t balance. The 
assumptions are not necessarily optimum for the range 
of vehicles considered, but they represent an average 
of the features of TRIESTE II, ALUMINAUT, and ALVIN 
as they apply to the family of vehicles considered. 

2, A computer program vas developed to optimize hull 
vei^t for a given depth. The program is described 
in Appendix A. Using the computed resiilts of hull 
vei^t versus depth as a base, the veight of payload, 
batteries, dropable ballast, main propulsion units, 
and "all other" variables was added. This resulted 
in a basic vehicle capable of proceeding to a given 
depth without the need of flotation material. At 
lesser depths additional ballast is required and at 
greater depths added lift is necessary. 
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3* Vehicles of various size were used in the con^jutation 
of hull weight versus depth. Size was varied in two 
ways: hy holding outside diameter constant while increas- 

ing length and hy holding the length to diameter ratio 
constant. The resxilts determine changes in size with an 
increase in depth for the two methods of increasing dis- 
placement. 

4, Finally^ flotation is considered as a means of hull 
support at increased depth. The addition of flotation 
necessitates an addition of the variables mentioned 
above^ batteries^ ballast^ etc. This reduces the effect 
of the flotation material. Curves were constructed for 
syntactic foam^ probably the best flotation material now 
available^ and for materials of lower densities which 
might be developed in the not-too -distant future. These 
curves take into accoTuat the weight of additional equip- 
ment to arrive at a reduced flotation effectiveness. 
Entering the curves with the additional lift required by 
a basic vehicle at a given depths one determines the 
volume of flotation material required^ wei^t of flotation^ 
and the wei^t of the added hardware necessitated by the 
additional displacement. The result of primary interest 
is the displacement of flotation to be added to the dis- 
placement of the basic vehicle. Total displacement is then 
plotted against depth for conparlson of size with similar 
curves for the larger pressure hulls considered. 
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B. Assuirg^tions 



1. A famly of cylindrical vehicles with hemispherical 
end caps is considered. 

2. A basic vehicle with outside diameter of seven feet 
and l/D ratio of 2,0, where L is the length including 
the hemispherical ends^ is considered. Hulls with 
eight-foot and 10-foot diameters^ also with L/D ratios 
of 2.0; are considered for comparison as are vehicles 
of seven-foot diameter with L/D ratios of 2.5 and 3-0* 
ThuS; size is varied in two ways; hy holding L/D constant 
and hy holding D constant. 

3* Pour materials are considered for construction as 
follows ; 



Material 


Yield Strength 


Modulus of Elasticity 


Density 


Steel 


80;000 p.s.i. 


30 X 10^ p.s.i. 


490 Ib./ft.^ 


Steel 


150;000 p.s.i. 


30 X 10^ p.s.i. 


490 Ib./ft.^ 


Titanium 


120;000 p.s.i. 


6 

15 X 10 p.s.i. 


281 Ib./ft.^ 


Alviminum 


60;000 p.s.i. 


6 

10 X 10 p.s.i. 


168 Ib./ft.^ 



Glass was also included for coniparison in the computation 
of hxill weights; using the same theory applicable to metals. 
This theory may; but probably does not give realistic values 
for the wei^t of glass hulls. The properties used for glass 
are as follows: 

Glass 3QQ;000 p.s.l. 14 x 10^ p.s.l. l84 Ib./ft.^ 
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The values of hull 'wei^t for glass are tabulated 
in Tables I - III but are not plotted. 

Pover requirements are based on a sustained speed 
of five knots for a period of ten hours. The auxiliary 
load is assumed to be 10 kilovatts. 

5. Main propulsion power is assumed proportional to the 
product of displacement to the two -thirds power and the 
speed in knots cubed. Since the speed is fixed at five 
knots ^ power is proportional to displacement to the 
two -thirds power. 

6. Silver-zinc batteries with the following properties 
are assumed for the power source. 

Capacity 6.5 kilowatt hours/cubic foot 
23.0 pounds /kilowatt hour^^^ 

7- Payload including complement and instruments is held 
constant at 2_,000 lbs. for all vehicles. 

8. Minimum volume reqviirement for payload and equipment 
was assumed to be 3^6 cu.ft. which is the internal 
volume of the seven-foot diameter, L/D = 2 pressure hull. 



(1) Keays, K. , "Parametric Study of Two Man Search and Work Submarines," 
General Dynamics/ELectric Boat contract study, page 24. 
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9. Auxiliary buoyancy is assured initially to be 



syntactic foam vitb a density of 44 lbs. per cubic 
foot. Buoyant materials vitb lesser densities of 
0.50 (32.lUper cu.ft.) and O.30 (21.33 cu.ft.) 
are st-udied for conparison. 

10. The remaining veights vbich do not fit into the 
categories of hifLl^ batteries^ or payload vary vith 
a number of parameters. Some vary vith displacement, 
some vith depth, some vith speed and some vith external 
or internal surface area. It is not the purpose of 
this paper to analyze each of these vei^t groups, but 
estimates must be made in order to concentrate on the 
buoyancy requirements. It is assumed, therefore, that 
the folloving veight categories in addition to hull, 
battery and payload vill be considered. 

Main Propulsion units - function of displacement - 

yfo of displacement 

iTotation material - function of buoyancy requirements 

and depth 

Ballast - function of total displacement - 8^ of displacement 

All Other - - . - - ^ 



(2) Wenh, E., Jr., DeHart, R. C., Mandel, P., and Kissinger, R., Jr. 
"An Oceanographic Research Submarine of Aluminum for Operation t 
15,000 Feet", Transactions R.I.N.A.., Vol. 102, Ho. 4. Oct. I960. 
Based on estimates made for ALDMZKAIJT, page I5. 
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C. Basic Vehicle Weights 



1. A computer program vas developed based on the class 
notes from course 13.15^ M. I. T., by J. Harvey Evans. 

The program optimizes frame size and frame spacing to 
arrive at a least vei^t solution for a hTJll of a 
given size, constructed of a given material. The com- 
puter program is described in Appendix A. The results 
are tabialated in Tables I - V, 

2. Figure VI is a graph of battery weight versus total 
displacement. For vehicles above 30,000 lbs., battery 
weight can be treated as a constant plus a percentage 
of the displacement over 30,000 lbs. with an error less 
than one-tenth of one percent of total hull displacement 
at all times. For a vehicle to which flotation material 
is to be added, the additional battery wei^t equals 
1.4 per cent of the additional displacement. 

3« Ballast, propulsion \mlts, and "all other" also are 
functions of displacement and, like batteries, these 
weights may be treated as a constant which is part of 
the basic vehicle, plus a weight equal to a percentage 
of the displacement of the bvioyant material added. Added 
wei^t due to an increase in displacement amounts to 
20,4^ of the added displacement. 
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4. Figures II - V are plots of veight versus depth and 
take into account hull vei^t^ batteries, ballast, 
propulsion units, "all other", and payload. Weights 
plotted in Figures II - V are exclusive of any flotation 
material. Displacement is also plotted in Figures II - V 
for the differently configured vehicles. Comparison of 
displacement to veight yields the ballast or the added 
lift required for neutral buoyancy. 

D. Buoyancy 

1. The buoyancy of a vehicle is equal to the vei^t of 
salt vater displaced. Flotation material must be added 
to a vehicle vhich veighs more than it displaces; ballast 
must be loaded aboard a vehicle vhich veighs less than 
its displacement. 

2. The addition of flotation material necessitates an 
addition in the requirements for batteries, dropable 
ballast, main propulsion units, and "all other". Having 
linearized battery vei^t as a function of added dis- 
placement, all of the above items may be treated as a 
single veight addition proportional to the flotation 
material displacement. 

3. The lift produced by flotation material is also a 
percentage of the flotation material displacement. The 
percentage may be treated as the effectiveness of the 
material. Syntactic foam, at 44 lbs. per cubic foot 
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is 31-3^ effective. A material vith a density 
vith respect to sea vater of O.5 is 50/^ effective 
and one vith a density of 0.3 is 'JO'fo effective. 

4. The effectiveness of the flotation material is 
reduced hy the additional weight described in 
paragraph 2 of this section. The reduced effective- 
ness is the flotation material effectiveness less 
the percentage of added weight. For the assumptions 
made in this study, the reduced effectiveness of 
syntactic foam is IO.9 (31*3 - 20.^1-); of the 0.5 
density material it is 29.6; and of the 0.3 
density material it is 49.6. Since this reduced 
effectiveness is the true indication of lift contri- 
buted to a basic vehicle, it can be seen that 
reducing the density from O.687 to 0.5 Increases the 
effective lift almost threefold. Flotation material 
displacement is reduced to one-third of the initial 
requirement and added batteries, dropable ballast, 
etc., is only one-third as much. 

5. Figures VTI - IX are graphs of flotation material 
characteristics and weight additions. Entering along 
the ordinate of the graph with the lift required by a 
basic vehicle at a given depth, the curves provide the 
following information: 
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Ctunre 1 - Lift applied to the basic vehicle 

Carve 2 - Lift applied to additional batteries, 
ballast , etc . 

Carve 3 - Total lift provided by the flotation 
raaterial 

Curve 4 - Dry "weight of floatation material 
to be added 

Curve 5 - Salt water displacement of flotation 
material to be added 



These curves apply to any of the vehicles in the family 
of vehicles considered, provided that flotation material is 
to be added to provide additional lift. 
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III. PRESENTATION OF RESULTS 



Pressure Hull Sizes 
Figure I 








Displacement of Hulls 




Hull Diameter 


Jj .O.A. 


l/d 


Displacement 
(LIs. Salt Water 


7^ 


14' • 


2.0 


28,650 




1T.5 


2.5 


37; 300 




21.0 


3.0 


45; 950 


8' 


i6.o 


2.0 


42,900 


10' 


20.0 


2.0 


83; 500 



- 12 - 



Table I 



Computed Cylinder Weights L/D = 1.5 



Collapse Depth 
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Table II 



Coniputed ¥ei^t of 2 Hemispherical End Caps 



Collapse Depth 
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Table III 

Comprcted Pressure Hull Wei^ts l/D = 2.5 
Collapse Depth 



D 


3,^00 ft. 


10,200 ft. 


17,000 ft. 


34,000 ft. 


Steel, HY-80 


T' 




22 ,l 40 


35,700 


68,100 


8' 


11,850 


32,700 


52,700 


101,000 


10' 


22,800 


62,800 


101,500 


193,500 



Steel, HY-150 



7 ’ 




12,530 


20,000 


37,900 


8' 




18,470 


29,400 


56,000 


10 ' 




35,400 


56,500 


107,900 



Titanium 



7 ’ 




9,380 


l 4 , 46 o 


27,080 


8' 


5,180 


13,470 


21,300 


40,250 


10 ' 


9,800 


25,900 


41,100 


77,600 



Aluminum 60 



7 ' 


3,580 




15,600 


30,200 


8' 


5,270 


14,300 


23,100 


44,700 


10 ' 


10,080 


27,400 


44,400 


86,500 



■ Glass 



7 ' 




3 ,o 4 o 


4,630 


7,900 


8 ' 




4,600 


6,840 


11,730 


10 ' 




9,300 


13,600 


22,500 
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Table IV 



Con5)uted Pressure Hull Weights L/D = 2.0 



Collapse Depth 




- l6 - 



Table V 



Computed Pressixre Hvill Weights l/D = 3*0 



Collapse Depth 
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Table VI 



Weights Added to Bare Huills to Arrive 
at Basic Vehicle Weights 

(Figures II-V) 



Diameter 


T’ 


l/d 


2.0 


Battery Wt. 
(Figure VI ) 


3; 300 


Ballast, Etc. 
(19/^ X Displ. ) 


5; ^50 


Payload 


2,000 


Total 


10;T50 




8' 


l/d 


2.0 


Battery Wt. 


3; 600 


Ballast, Etc. 


8,150 


Payload 


2,000 


Total 


13,T50 



7' 


7' 


2.5 


3.0 


3;550 


3,750 


7,100 


8,7^0 


2,000 


2,000 


12,650 


lk,k^0 



10 ' 

2.0 

i^,300 

15; 900 

2,000 

22,200 
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Figure XII 



Sample Computer Solution 

Material - Alumimmn 60 
Collapse Depth - 17_,000 feet 
Hull Diameter - 7* 
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rv. COKCLUSIOKS MD ESCOM-IEKDATTOHS 



A. Figure X displays the results required by section I of this thesis. 
Overall vehicle size is plotted against depth for hulls constructed of 
four different materials. The effect of increasing depth by changing 
hull size or by adding flotation material with pressure hull size con- 
stant is cong)ared •vd.th the following results. Increasing hull size by 
varying vehicle length provides greater depth than building a larger 
vehicle with L/D constant for HY-80 and Titanium. For aluminum the 
opposite is true and for HT-150 both methods have eq.ual effect. Therefore^ 
the results with respect to a preferred method of increasing size are 
inconclusive^ one method having nearly the same effect as the other. 
However, the addition of flotation material is to be preferred over 
changing vehicle size by either method for all but HY-80, assuunlng 
syntactic foam is the material used. For depth Increases greater than 
about 5; 000 feet, a smaller vehicle results if flotation material is 
added. And, of course, the cost of a given volume of flotation material 
will be much less than the cost of increasing hull size by an equal 
amount. 

Recommendation: Deep submergence pressure hulls should be 

built large enough to satisfy minimum volume requirements. If added 
lift is necessary, it should be provided by flotation material. 
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3. T'he effect of improving the characteristics of flotation materials 
is graphically demonstrated in Figure X. If a flotation material could 
he produced vith a density of O. 5 O, there vould he no question about the 
preferred method of increasing depth. 

Recommendation: As flotation material of lower density is 

developed^ its use should he considered in the desiga of deep suhmersihles , 
Costs will he cut and smaller^ less complex deep suhmersihles will resTilt. 

C. Although the results concerning a preferred method of changing vehicle 
size were inconclusive^ it can he said that there is prohahly an optimum 
hull configuration. With the computer program developed, one could arrive 
at hifLl proportions for a least weight vehicle with a given displacement. 

At first glance it might appear that this would always result in a sphere, 
hut with conventional framilng or sandwich hull construction, this prohahly 
will not he the case. The program can also he used to deteriaine optimum 
frame configuration, and it appears to he a very promising tool in evalu- 
ating current desiga practices as they apply to vehicles designed for 
greater depths. 

Recommendation; The computer program described in Appendix A 
shoifLd he used to continue the study of hull configuration started in this 
paper and it coifLd he used as well for a study of the parameters of sub- 
marine design and their interaction on one another. 

D. Figure XI is a graphic illustration of the results obtained by using 
the methods outlined in the notes from course 13-15 an<i from earlier cal- 
cvilations based on unstiffened cylinder construction. Some difference 
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exists betveen the wei^t of unstiffened cylinders and the weight of 
the stiffened cylinder as predicted by the program in Appendix A, 

It is natural that a difference exists but in some cases the unstiff- 
ened cylinder is limited to a depth 38 ^ less than that of a stiffened 
cylinder of eq,ual weight. This large difference is difficult to explain 
and it may be an indication that the notes on which the computer program 
is based may give incorrect collapse depths when considering deep sub- 
mersibles^ or_, that the con^juter program contains some flaw. However^ 
the program was tested using as inputs data from relatively shallow- 
depth submersibles^ and it duplicated almost exactly the calculations 
which resulted when those vessels were designed. 

Recommendation': The computer program^ probably the most 

valuable by-product of this thesis^ should first be checked once again 
to ensure that the formulae have been translated into computer language 
properly^ and then the output should be compared with experimental results 
to determine: 

1, the validity of the theory with respect to 
conventional steel construction, 

2, the applicability to other metals, and 

3 , the application of the theory to vessels 
with t/D ratios approaching 0.1. 

Having checked the program and its general applicability, it will be 
possible to determine the necessity of modifying the basic theory for 
metals other than steel or, more probably, for deep submergence vehicles 
with very thick shells. 
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APPEI'IDIX A 



Computer Design Program to Optimize Submarine Hull Wei^t^ based on 
class notes of Course 13.15^ M. I. T.^ by J. Harvey Evans ^ Professor 
of ITaval Architecture. 

1 . Background 

a. Failure of a stiffened shell occurs in one of three 
modes. The first occurs if the shell buckles between 
the frames while the frames retain their circularity. 

This is referred to as elastic or plastic instability. 

The second mode of failure inrolYes yielding of the 
shell. In this mode the shell remains circular but 
the yield strength of the material is exceeded. The 
third mode of failTore occurs when the shell and frames 
fail together and this is known as general instability. 
Ideally, for optimizing strength with least weight, a 
shell should fail in all three modes simultaneously. 

This program optimizes three parameters, shell thickness, 
frame area, and frame spacing in an attempt to find the 
least wei^t solution. 

b. In such a solution it would be ideal to optimize the 
dimensions of the frame as well as the other parameters 
mentioned, but this would have introduced another degree 
of complexity to an already complex analysis. Instead, 
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i 





a frame typical of those presently \ised in submarine 
construction uas selected as a parent and all frames 
are scaled from the parent member. The cross section 
of the parent member is shovn in Pigure A-1 below. 

All dimensions are with respect to veb 

thiclmess. 




w 



c. Having selected a parent frame member^ it is only 
necessary to settle on the length of cylinder between 
bulkheads and the mean cylinder diameter in order to 
applythe formulae presented in the set of notes refer- 
enced above. The general method of approach outlined 
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is that of selecting a set of scantlings and checldLng 
their applicability. An understanding of the paper is 
necessary for a full understanding of the program and 
its capabilities. However, lack of such understanding 
will not prevent successful use of the program. 

. The program is based on the solution of six eqviations. 
These are; 

1. Hoop stress formula which predicts a pressure 
at which shell yield occurs. 

Per = 2(t/c) CT 

y 

This formula is used to determine shell thiclmess 
for a contecplated design. Shell thickness, 

. Per X D 

Experiments indicate that the predicted thiclcness 
may be reduced in some cases by as much as 20^. 

The computer program, however, corrpletes two 
solutions , at predicted ’t ' and at 95^ of 't ' but 
a minor modification to the program allows 't' to 
be reduced to any given percentage of the value 
predicted by the hoop stress formula, or it allows 



ji 

J 




'f to be reduced in increments until a least 
vei^t solution is found wMcb satisfies all 
other criteria but ignores the hoop stress 
prediction. 

2. Pormula 92 a modified predicts plastic insta- 
bility failure. 



D 



m 



1 



' cr 



1 + H' 



0.85 - B 1 + o.l48(\ - s)'" 



1 + 3 



The critical pressure it predicts must be equal 
to or greater than collapse depth pressure, (see 
page h2 for a list of symbols and their meanings.) 

3. ¥enk’s formula predicts collapse pressvire based 
on the circumferential stress at the shell mid- 
thichness midway between frames. 



P 

cr 



2ta 

z 



D 



X 



1 




4. Lunchick predicts "yield pressure" based on a 
three hinge failure mechanism. 



P 

y 



2tcr 

3 

D 



X 



"Vx^ - xz + 



This formula is normally within 3 - 55^ of the collapse 
pressure predicted by ¥enk. 
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5» Bryant has taken a collapse prediction by 
Kendrick and put it into a useful form. This 
forimfLa is based on the number of conp>lete cir- 
cumferential waves formed by the plating and 
frames buckling as a unit. It is solved^ normally 
for 2., 3, or 4 waves and the first critical mode 
of failure (i.e., the one for which the critical 
load is smallest) is the value sought. 



T3 2Et 

P = X 

cr D 

m 



m 



/2 m _\/2 ^ 

(n + -jp - l)(n + m ) 

Formulae 2^3/4 and three variations of 5/ predict 

critical press-ures. It is necessary to solve these 

formulae and then select that answer which is lowest 

as the critical one in the case under consideration. 

6. The last formula used predicts the stress in the 
frame flange due to the combined effects of total 
load including shear and bending due to out of 
roundness. This formula is based on the critical 
pressures predicted by Bryant’s formula j thus^ it 
is also solved three times in the program. 



8 . (n - 1) 

2.2 



max 




1 + 
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The highest values predicted hy this formula 
must not exceed the yield stren^h of the material. 
Thvis^ the hipest value is selected hy the program 
as the critical value. 

e. The program solves for shell thiclmess -using the hoop 
stress formula and then computes six solutions for critical 
pressure. The lovest of the six is termed the critical 
pressure for solution of the problem. Then three solutions 
of maxinrum stress are made^ and the highest of these is 
selected as the critical stress. At this point the program 
computes the amount hy which critical pressure must he 
changed -to agree -with co2JLapse pressure, and the amiount hy 
which the highest stress must he altered so it does not 
exceed yield strength of the material. If the program can 
he made to alter frame size and frame spacing until these 
■values agree; critical press-ure with collapse press-ure, and 
maximum stress with yield stress, then for a given shell 
thickness, a least weight solution will he possible, 

f , The solution is fo-und by solving the problem with frame 
area incremented a \uilt amo-unt and with frame spacing incre- 
mented a unit amount. These two solutions yield changes in 
the critical values of pressure and stress; that is, they 
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give the partial differentials of pressure and stress 
vith changes in frame area and frame spacing. It is 
then a sinpile matter to assume that linearization is 
possible and solve two equations in two unlmowns for 
the number of unit changes to be made in frame area^ A_, 
and frame spacing^ SP. 

dA + ^ dSP = (P -P .,) 

gA oSP ^ collapse crit*^ 

II dA + -41 dSP = (S . , , - S .,) 
oA oSP ' yield crit"^ 

Solving for the two unlmowns dA and dSP_, the original 
estimate of frame area and frame spacing can be changed. 
The problem is solved again with the new values of 
flame area and frame spacing and if the critical pressure 
is within a given percentage (the program specifies k^o 
on the safe side) and if the critical stress is also safe 
and within a set tolerance (these tolerances can be 
changed if a higher degree of accuracy is desired),, the 
program will cong>ute weight of compartment including shell 
and frames. (One word of warning. In order to ensure a 
satisfactory solution in all cases,, it was necessary to 
make frame spacing continuous. Therefore^ the nxunber of 
frames, eq,ual to compartment length over frame spacing 
may not be a whole number. It will be necessary to select 
the next higher nixmber of frames to ensure a safe solution 
with a subtly hi^er wel^t than that computed.) 
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g. If, at the end of the first iteration^ the values of 
predicted collapse pressure and predicted stress are not 
vithin tolerance^ the program vill continue to cycle 
through the series of calculations inproving the solution 
until it is vithin tolerance. The program then increments 
shell thickness to 95^ of its predicted value and solves 
the problem again. The veight determined is compared vith 
the vei^t determined vith the original value of 't* and 
the "least veight" is printed before the problem goes on 
to the next set of data. 
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LIST OF SYMBOLS 



Some of the symbols in the folloving list refer to variables 
vhich are used to simplify the basic formulae. When such variables do 
not have explicit definitions, reference is made to a formula number, 
figure number, or page number, vhere the symbol is used in the notes 
by Professor Evans. 



A 


Cross sectional frame area 




Area of one frame plus shell under the veb 


b 


Web thickness of the ring frame 


B 


Page l8 


c 


Formula 31 


D 


Outside diameter of hull 


D 

m 


Diameter of hull to mid-shell 


E 


Young's Modulus 


F 


Page 31 


H 


Figure 8 


^1 


Moment of inertia of a T section 




Frame Spacing 


m 


3.1^ D 1 2L (rounded off to the nearest vhole number) 

• 


n 


2, 3, or 4 - no. of lobes for the three modes of 
failure by general instability 


P ■ 


Predicted collapse pressure 


^cr 


Pressxore at vhich yield occurs 



k2 - 



Lovest collapse pressure of the 6 pressures calculated 



^crit 

P Yield press-ore - same as P 

y cr 

P^ Pressure at collapse depth 

R Radius of hull to mid-shell 

m 

^crit stress of the three computed 

SP Frame spacing 

t Shell thickness 

X Page 3 ^ 

z Page 34 

P Page 18 

e Page 9 

A Slenderness ratio - Figure 4 
YS Yield stress 



2 . Input 

a. Input consists of the length and mean diameter of the 
cylinder under consideration, the characteristics of the 
miaterial to he used in construction, including yield 
strength. Young's modulus and density of the material, 
the collapse depth, allovahle out-of-roxmdness of the 
hull, and the hest estimate available of frame area and 



frame spacing. 






Inputs 


Units 


Definition 


A 


. 2 

in 


Area of frame 


CD 


feet 


Collapse depth 


CL 


inches 


Compartment length 


D 


inches 


Mean diameter of hull 


Delta 


Inches 


Allovable out -of -roundness 


Dens 


Ib/ln^ 


Density of material 


E 


p.s.i. 


Young's modulus 


SP 


inches 


Frame spacing 


SY 


p.s.i. 


Yield strength 


(Delta, allovable 


out -of - 


-roundness of the hull, should 



normally be limited to one-half inch or half the thick- 
ness of the shell. For heavy shell plating this may be 
increased, but for shell plating less than one-inch 
thick Delta vill be altered by the program to one-half 
the thickness of the shell plating. ) 
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b. All input for a single problem is placed on one card 
in the following form: 



A 



7 



SY 



17 



E 



27 



33 



CD 



39 



CL 



^5. 



Dens ^^SD ^^Delta 



The numbers indicate the card space number at which each item 
should start. Sample input follows; 



70 60000. 30000000.-360.— 400.-540. — .2875-30. -.50- 



The sanpile input indicates the following: 



Frame area 


70 sq.in. 


Yield Strength 


60,000 p.s.i. 


Young's Modulus 


30 X 10 p.s.i. 


Mean diameter of hull 


360 inches 


Collapse depth 


4^000 feet 


Compartment length 


540 inches 


Density of material 


.2875 It** in 


Frame spacing 


30 inches 


Allowable out-of -roundness 


0.50 inches 
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An input deck, ready for the computer "will consist of the 
follo'wing cards : 



Meaning 


Space 


Emher 




1 


T 


Identity 


* 


M3999-3504,EMS;DEBUG, l, 500, T5 Peterson 


Execute 


* 


XEQ 


Binary 


* 


BINARy (or PORTRAE) 


Comment 


C 


Peterson (this card included only if a 

fortran deck is used) 






First binary or fortran card of program 



Last card of program 





El© 


* 


DATA 




First data card 
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m 



3 



MINIMUM HULL WEIGHT PROGRAM 



DIMENSION PCAM (6),PCW(6)»PCL(6)*PCB1(6) ,PCu2 ( 6 ) »PCB3 ( 6 
1 ) ( 6 ) »SMK2 (6 ) »SMK3( 6 )»P(10)»S(4),W(40).PDIFS(10) 

2,SDIFS( 10) 

10 READ 1340»A»SY»E,D.CD»CL»DENS»SP*DELTA 
PCD = CD^-64. /144. 

PRINT 1010, SY 
PRINT 1040, CD 
PRINT 1090, PCD 
T = CD*64.*D*0. 5 / ( 144.->SY ) 

IF(T-0.25)20,30,30 
20 T=0.25 

30 IF{T-1.0)40»50,50 

40 DELTA=T/2. 

50 continue 

SY=0,98*SY 

coai.oa-^CD 

PRINT 1020, E 
PRINT 1030, D 
PRINT 1050, CL 
PRINT 1060, DENS 
PRINT 1080 , DELTA 
PCD=1.02^PCD 
LR = 1 

15 L = 1 

L0=1 
LP = 1 
NP = 1 
NQ=1 

C PRELIMINARY CALCULATIONS 

60 FN=CL/SP+1. 

IF(A) 64,64,65 

64 PRINT 1390, A 
GO TO 10 

65 IF(SP)66.66,67 

66 PRINT 1400, SP 
GO TO 10 

67 LO=l 

TW=SQRTF( A/19. 18) 

THETA = 18.2-*^SP/ (D*SQRTF( lOO.^T/D) ) 

SINH1 = ( 2.71828**THETA-1./2.71828*^;-THETA) /2.0 
C0SH1 = ( 2.71828**THETA+l./2 .71828^i-*THETA /2. 0 
GK=(SINH1-S INF (THETA) ) / ( 8 1 NHl + S I NF ( THETA ) ) 
GN={CGSH1-C0SF( THETA) ) / ( S I NHl+S I NF ( T hETA ) ) 
THA2=THETA/2, 

SI NH2=( 2.71828**THA2-1 ./2, 71828**THA2 ) /2 . 

COSH 2= { 2.7182 8^<-*THA2 + 1./2.71828**THA2 ) / 2. 

GH = -(3,*SINH2^^COSF(THA2)+COSH2*SINF(THA2) ) / (SINHl + SINF 
1 { THETA) ) 

Q=-( SINH2*COSF (THA2 )+C0SH2*SINF( THA2) ) / ( S I NHl+S I NF ( T hE 
ITA) ) 

GT=-( SIN,-I2*C0SF ( THA2 ) -COSH2*S I NF ( THA2 ) ) / ( S I NHl + S I NF ( TH 
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lETA) ) 

30A= ( ( SP/D ) ^-> 1 / ( T/D ) ^<-*3 ) *^-0.2 5-K-SQRTF( SY/E ) 
EPSL=(0»A5^<-SY'«-D)/(E*T ) 

3 = TW*T/ ( A + TV.'*T ) 

BETA= ( 1 1 .•!<-GN^^T^-*2 ) / ( SQRTF( 100 . ^<-T / D ) * ( A+TW*T ) ) 

BM = 3 . 1A16*D/ ( 2.«-CL ) 

IBM=XFIXF ( BM+0.5 ) 

BM=FLOATF( IBM) 

F=1 ,+2,*Q* ( 0.85-B ) / ( 1 .+BETA ) 

X=1 ,+GH* ( 0. 85-B ) / ( 1 ,+BETA ) 

2 = 0. 5 + 3.6 2^>‘GT*( 0.8 5-B) /( 1.0+BETA) 

EL = 1 .57*SORTF( D*T/2. ) 

IF (EL-SP) 80 »80 >70 
70 EL=SP 

8 0 ERT IA1 = 102.2*TW**4*( 3.*EL*T / ( EL^^ T + 1 9 . 1 8*TW**2 ) * ( 0 . 1 25* 
lT/TW+1.43)**2+0.79) 

AB=A+TW*T 

C= (43 . 5*TW**3+ ( T*EL ) * ( 8 .0*TW+0.5*T ) ) / ( 19 . 18*TW**2+EL*T ) 
PCAM(L)=2.*T*SY/(0*( l.+GH*( ( 0.85-B) /( 1.+8ETA) ) )*( l.+O. 
1148*( ( 80A**2-EPSL) **2 ) ) ) 

N0 = 1 

P( NO ) =PCAM ( L ) 

NO=NO+1 

PCW ( L) =2.*T*SY/ (D*SQRTF ( F**2-F/2 . +0 . 2 5 ) ) 

P(NO)=PCW( L) 

N0=N0+1 

PCL ( L) = 2,0*T*SY/ ( D*SQRTF ( X**2-X*Z+Z*2 ) ) 

P( NO) =PC^( L ) 

N0=N0+1 
BN = 2 . 

PCB 1 ( L ) =2. *E*T*BM**4/ ( D* (BN**2+ ( . 5*BM ) **2-1 . ) * ( BN**2+B 
1M**2 ) **2 )+8.*( BN**2-1 . )*E*ERT lAl/ ( D**3*SP ) 

P (NO ) =PCB1 ( L ) 

N0=N0+1 
BN = 3. 

PCB2 ( L ) =2.*E*T*BM**4/ ( D* (BN**2+( .5*BM ) **2-1 . ) *( BN**2+B 
1M**2)**2)+8.*(BN**2-1. )*E*ERTIA1/ (D**3*5P) 

P ( NO) =PCB2 ( L ) 

N0=N0+1 
BN = 4. 

PCB3 ( L ) =2 .*E*T*BM**4/ ( D* (BN**2+( . 5*BM ) **2-1. ) *( BN**2+b 
1M**2 )**2 ) +8.*( BN**2-1 . ) *E*ERTIA1/ ( D**3*SP ) 

P(N0)=PCB3(L) 

BN = 2. 

SMKl (L) =(PCD*(0.5*D+0.5*T)*TW/ (AB*( 1.0+BETA) ) )*(1.0+(0 
1 .85*BETA)'/B ) 

2+E*C*(BN**2-1.0)*DELTA*PCD/( ( 0 . 5*D ) **2* ( PCB 1 (L)-PCD) ) 

N0 = 1 

S ( NO ) =SMK.l ( L ) 

N0=N0+1 
BN = 3. 

SMK2 ( L) = ( PCD*( 0.5*D+0.5*T )*TW/ (AB*( 1.0+BETA) ) )*( 1.0+(0 
1 .85*BETA) /B ) 

2 + E*C* (BN**2-1 . o') *DELTA*PCD/ ( ( 0 . 5*D ) **2* ( PCBl ( L) -PCD) ) 
S(NO ) =SMK2 ( L ) 
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c 



c 



90 

110 



100 

120 

122 

125 

130 

140 



150 

160 

162 

163 

164 



165 

170 



180 



190 



NO=NO+l 
3N = 4. 

SMK3 ( L ) = ( PCD*( 0.5^:-D + 0. 5*T ) *T W / ( AB ^ ( 1. 0 + 3 ETA ) ) ) * ( 1 .0+ ( 0 
1 .85*BETA) /B ) 

2 + E*C* ( BN** 2-1 . 0) *DELTA*?CD/ ( ( 0 . 5*D ) **2* ( PCSl ( L) -PCID ) ) 

S (NO ) =SMK3 ( L ) 

COMPUTE VOLUME OF MATERIAL IN ONE COMPARTMENT LENGTH 
V=0.78 54*CL*( (D+T/2. )**2-(D-T/2. ) *->-2 ) +0 . 78 5 4* ( ( D-T/2* ) 
1**2- ( D-T/2.- 

2 5.9 3*TW) **2 )*TW*FN + 0.78 54*-FN*6,4*TW*( ( D-T/2 . -5 . 93*TW ) * 
3*2- ( D-T/2,-8 .0*TW) **2 ) 

COMPUTE WEIGHT OF MATERIAL IN ONE COMPARTMENT LENGTH 
W(L)=V*DENS 
L = L + 1 

GO TO (90,165) ,LP 
DO 110 N=l,6 
PDIFS(N)=PCD-P(N) 

CONTINUE 
PDIFCR=PDIFS ( 1 ) 

DO 120 N=2,6 

IF (PDIFS(N)-PDIFCR) 120,100,100 
PDIFCR=PDIFS (N ) 

NP = N 

CONTINUE 

IF( 0.02*PCD-ABSF(PDIFCR) )130,125,125 

LQ=LO+l 

DO 140 N=l,3 

SDIFS(N) =SY-S(N) 

CONTINUE 
SDIFCR=SDIFS ( 1 ) 

DO 160 N=2,3 

IF ( SDIFS ( N ) -SDIFCR ) 150 ,160 , 160 
SDIFCR=SDIFS (N ) 

NQ = N 

CONTINUE 

IF ( 0.02*SY-ABSF (SDIFCR) ) 164,163,163 
LO=LQ+l 

GO TO ( 164,164,200) ,LQ 
LP = 2 

PCRIT=PCj-PDIFCR 
SCRI T=SY-SDIFCR 
GO TO ( 170,180,190) , LO 
AT = A 

A=0.96*A 
LO = 2 

GO TO 60 

DELPE1=PCRIT-P(NP) 

DELS1=SCRIT-S(NQ) 

A = AT 
SPT=SP 
SP=1 . 04*SP 
LO = 3 

GO TO 60 

DELPE2=PCRIT-P (NP) 

DELS2=SCRIT-S(NQ) 
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t 



200 



SP=SPT 

D0T=DELS2-DELPE2*DELS1/DELPE1 

D0P=SDIFCR-PDIFCR*DELS1/DELPE1 

PRY=DOP/DOT 

PRX= ( PDIFCR-PRY*DELPE2 ) /DEL PEI 

A=(1.+0.04*PRX)*A 

SP= ( 1 .-.OA^^PRY )*SP 

GO TO 15 

PRINT 1000, A 

PRINT 1070, FN 

PRINT 1140,PCAM(L-1 ) 

PRINT 1 1"0 ,PCW (L-1 ) 

PRINT lloO,PCL (L-1 ) 

PRINT 1170 ,PCB1 (L-1 ) 

PRINT 1200, PCB2 (L-1 ) 

PRINT 1210,PCB3(L-1 ) 

PRINT 1220 ,SMK1 (L-1 ) 

PRINT 1230,SMK2(L-1 ) 

PRINT 1240 ,SMK3(L-1 ) 

PRINT 1100, SP 
PRINT 1110, T 
PRINT 1120, TW 
PRINT 1180, V 
PRINT 1190,W(L-1) 

PRINT 450, THETA 
PRINT 450,SINH1 
PRINT 450,C0SH1 
PRINT 450, GK 
PRINT 450, GN 
PRINT 450,THA2 
PRINT 450,SINH2 
PRINT 450,C0SH2 
PRINT 450, GH 
PRINT 450,0 
PRINT 450, GT 
PRINT 450>BDA 
PRINT 450,EPSL 
PRINT 450, B 
PRINT 450, BETA 
PRINT 450, BM 
PRINT 450, F 
PRINT 450, X 
PRINT 450, Z 
PRINT 450, EL 
PRINT 450,ERTIA1 
PRINT 450, AB 
PRINT 450, C 
PRINT 1130, PCRIT 
PRINT 1130,SCRIT 
PRINT 1130,PDIFCR 
PRINT 1130, SDIFCR 
PRINT lUO, DELPEl 
PRINT 1130, DELPE2 
PRINT 1130, OELSl 
PRINT 1130, 0ELS2 
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t 



210 



PRINT 1130, DOT 
PRINT 1130, OOP 
PRINT 1130, PRY 
PRINT 1130, PRX 
GO TO ( 210,220,250 ) , LR 
WLEAST=W( L-1 ) 

PRINT 1370, WLEAST 
LR = 2 
TR = T 

T=0.95*T 
GO TO 15 

220 IF(WLEAST-W(L-1 ) ) 240, 240*2 30 
230 WLEAST=W( L-1 ) 

240 PRINT 1370, WLEAST 
T=1.04*TR 
LR = 3 

GO TO 15 

250 IF(WLEAST-W(L-1 ) ) 2 70 , 2 70 > 260 
260 WLEAST=W(L-1 ) 

PRINT 1370, WLEAST 

T=1.04^^T 

GO TO 15 

270 PRINT 1370, WLEAST 
GO TO 10 

450 FORMAT(F20.5) 

1000 FORMAT ( IHl , lOHFRAME AREA » 1 9X , 1H= , F i 0 . 2 , 3X , 7 hSQ . IN.) 
1010 FORMAT! 1H1,14HYIELD STRENGTH , 1 5X , 1H= , FI 0 . 2 , 3X , 6HP . S . I . ) 
1020 FORMAT! IX, 14HY0UN6S MODULUS , 1 2X , 1H= , F 1 3 . 2 , 3X , 6HP . S . I . ) 
1030 FORMAT! IX, 21HMEAN DIAMETER OF HUL L , 8X , 1H= , F 1 0 . 2 , 3 X , 6H I 

NCHES) 

1040 FORMAT! IX, 14HC0LLAPSE DEPTH , 1 5X , 1H=- , F 10 . 2 , 3X , 4HFEET ) 
1050 FORMAT! IX, 16HBULKHEAD SPAC I NG , 1 3X , 1 H= , F 1 0 . 2 , 3X , 6H I NCHE 

S) 

1070 FORMAT ! IX, 30HNO. FRAMES BETWEEN BULKHEADS =,F10.2,3X,6 

Hi-RAMES) 

1060 FORMAT I IX, 19HDENSITY OF MAT ER I AL , lOX , 1H= , F 1 0 . 4 , 3X , 1 6HL 

3. PER CU. I 

INCH) 

1080 FORMAT ! IX, 30HALLOWABLE OUT-OF-ROUNDNESS =,F10.2,3X,6 

HI NCHES ) 

1090 FORMAT ! IX, 30HPRESSURE AT COLLAPSE DEPTH =,F10.2,3X,6 

HP.S. I . ) 

1100 FORMAT! IX, 13HFRAME SPAC I Ng , 16X , 1H= , FI 0 . 2 , 3X , 6H I NCHES ) 
1110 FORMAT ! IX , 15HSHELL T H I CKNESS , 14X , 1 H= , FI 0 . 2 , 3X , 6H I NCHES ) 
1120 F0RMAT(1X,13HWEB TH I CKNESS , 16X , 1H= , Fl 0 . 2 , 3X , 6H I NCHES ) 
1130 FORMAT ! F20. 5 ) 

1140 FORMAT ! IHO ,30HCOLLAPSE7PRESS. FROM 92AMOD. =,F10.2,3X, 

6HP.S.I.) 

1150 FORMAT! IX, 30HCOLLAPSE PRESS. FROM WENK =,F10.2,3X,6 

HP.S. I . ) 

1160 FORMAT! IX, 30HCOLLAPSE PRESS FROM LUNCHICK =,F10.2,3X,6 

HP.S. I . ) 

1170 FORMAT ! IX, 30HCOLL PRESS FR BRYANT, 2 LOBES= , F 10 . 2 , 3X , 6 

HP.S. I . ) 

1180 FORMAT! IX, 30HVOLUME OF ONE COMPT MATERIAL =,F10.2,3X,1 
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2HCUBIC INCH) 



lES) 

1190 FORMAT! 1X,21HWEIGHT OF cOMP AR TMENT , 8X . 1 H= . FI 0 . 2 , : X » 6HP 

OUNDS ) 

1200 FORMAT ( lx .30HC0LL PRESS FR BRYANT, 3 LOB£S= , F 10 . 2 , 3X , 6 

HP.S. I . ) 

1210 FORMAT! IX, 30HCOLL PRESS FR BRYANT, 4 L03ES= , F 10 . 2 , 3X , 6 

HP.S. I . ) 

1220 FORMAT! IX, 30HMAX STR FR BUSHIPS, 2 LOBES =,F10.2,3X,6 

HP.S. I . ) 

1230 FORMAT ! IX ,30HMAX STR FR B'JSHIPS, 3 LOSES =,F10.2,3X,6 

HP.S. I . ) 

1240 FORMAT ! IX, 30HMAX STR FR BU'HIPS, 4 LOBES =,F10.2,3X,6 

HP .S. I . ) 

1340 FORMAT !F6.2,2F10.0,3F6.0*F6.4,F4.0,F4.2) 

1370 FORMAT! IX, 30HCOMPARTMENT OF LEAST WEIGHT =,F10.2,3X,6 

HPOUNDS) 

1390 FORMAT! IX, 30HACTION TERMINATED, NEGATIVE A= , F 10 . 2 , 3X , 7 

HSQ. IN.) 

1400 FORMAT! IX, 30HACTION TERMINATED, NEGAT. SP =,F10.2,3X,6 

HINCHES) 



END 



NOTE 



FORMAT STATEMENTS REQUIRE A CARD FORM 72 CHARACTERS 
IN WIDTH. IN ORDER TO FACILITATE A NARROWER PRINTOUT 
THOSE FORMAT STATEMENTS LONGER THAN 60 CHARACTERS WERE 
BROKEN AND THE EXTRA CHARACTERS PRINTED ON THE FOLLOWING 
LINE WITH AN EXAGERATED INDENTATION. A USER COPYING THE 
PROGRAM SHOULD TYPE, FOR EXAMPLE, LINE 1400 AND THE NEXT 
INDENTED LINE ON THE SAME CARD. 
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I 



I 




I 





p 



COI-IPOTER EROGRAI-1 SOBOLS 



A 


Cross-sectional area of a frame 


AB 


Cross-sectional area of a frame pl\is the area of 
shell voider the veh 


AT 


Same as A 


B 


Page l8 


EDA 


Lambda - Figure k 


BETA 


Page l8 


EM 


m (see formula symbols) 


BN 


n (see formula symbols) 


C 


Formula 31 


CD 


Collapse depth 


CL 


Compartment Length 


COSHl 


COSH e 


C0SH2 


COSH e/2 


D 


Hull diameter 


DELPEl 


P . . - P . . (after A is incremented) 


DELPE2 


P . , - P . , (after SP is incremented) 


DELSl 


^crit ~ ^crit ^ incremented) 


DELS2 


^crit ~ ^crit incremented) 


DELTA 


Allowable out-of -roundness of the hull 


DENS 


Density of construction material 


DOT, DOP 


Intermediate values in the solution of two 
differential equations 


E 


Young ^s modulus 
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EL 


Effective length of shell plating 


E?SL 


c, page 9 


ERTIAl 


Koment of inertia of a T section 


E 


Page 31 


FN 


ilo. of frames in one compartment length 


ffil 


Fig\jre 8 


GK 


Figin’e 8 


GN 


Figure 8 


GT 


Page 34 


P 


Predicted collapse pressvire = PCAM_, PCW, PCL, PCBl, PCB2, 
PCB3 


PCAM 


Predicted collapse pressure from modified formula 9^A 


PCBl 


Predicted collapse pressure from Bryant's (Kendrick) 
formula (1st failirre mode) 


PCE2 


(2nd failure mode) 


PCB3 


(3rd failiare mode) 


PCD 


Pressure at collapse depth 


PCL 


Predicted collapse pressure from Lunchick 


PCBIT 


Lowest collapse pressure of the 6 predicted 


PCW 


Predicted collapse press-ure from Wenk 


PDIFCR 


Highest of the 6 PDIFS values 


PDIFS 


= PCD - P (six of these calculated) 


PRX 


No. of increments hy which SP must be changed in 
reaching an optimum solution 


PRY 


No. of increments by which A must be changed in reaching 
an optimum solution 


Q 


Page 31 
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1 



I 







k 



n 



f 



s 

SCRIT 

SBIRCR 

SDIFS 

SIRHl 

SI1JK2 

SIW. 

SMK2 

SMK3 

SP 

SPT 

SY 

T 

THA2 

THETA 

TW 

V 

W 

WLEAST 

X 

Z 



Predicted stress = SMECL^ SMK2^ SMK3 
Highest stress of the 3 predicted 
Lovest of the 3 SDEFS values 
SY - S (three calculated) 

sum © 

SINH e/2 

Predicted stress (2 lobe mode of fail\ire) 

Predicted stress (3 lobe mode of failure) 

Predicted stress (4 lobe mode of failure 
Frame spacing 
Same as SP 
Yield stress 
Shell Thickness 
0/2 

0 (Figure 8) 

Web thickness 

Volume of material in shell and frames over one 
compartment length 

Weight of material in shell and frames over one 
compartment length 

Weight of the lightest compartment of those considered 
Page 34 
Page 34 
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b 




4 . Output 



a. Given the folloving Inputs for a section of pressure 
hull: 



1 .) 


Diameter 


2 .) 


Length between transverse bulldaeads 


3.) 


Collapse depth 


^.) 


Construction mterial 


5.) 


Tolerance to which the construction 
material will be fabricated 


the computer 


program produces recommended values for 


1 .) 


Shell thickness 


2 .) 


Erame size, and 


3.) 


Erame spacing. 



Wiihthe above information, it is possible to compute 
displacement of the cylinder and weight of material used 
in construction of the hull. The coirputer program does, 
in fact, calculate weight of hull since it was decided 
that "least weight" would be the criterion for optimum 
design. 

b. But the process of design involves more than one answer 
to one problem. Normally, a designer will study several 
vehicles with various characteristics and he will attempt 
to gather. as much information as possible about each one. 
It is not enough to know that a computer program has 
recommended certain values for a certain vehicle. Each 
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M 




'one recoamended vs.lu.es miisi; le checked for ac curacy 
logic ^ and applicability to the planned use and the 
construction methods currently available. There is 
alvays the possibility that a coD^iuter program vill 
recommend as optimum a set of values impossible to adhere 
to in actual construction. 

c. To aid a designer in comparing vehicles and to aid in 
analysing tha computer programi Itself, or possibly to 
help in studying the design parameters of submarine con- 
struction, the computer program listed in this section 
prints out a large number of values computed in the 
preliminary calculations as veil as the final ansvers 
to the main problem. Once a user has gained familiarity 
vith the program he vill probably eliminate much of the 
printed output by vithdraving from the input deck those 
PREST statements referring to the unnecessary values. 

He may also add some cards or some additional calculations 
vhich vill help him in his design or analysis. 

Referring to the sample output (page ^0 the folloving 
values are printed out after each solution at a . given 
thickness. 
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Computer Output 
Page One 



1. Yield strength of the material 

2 . Collapse depth 

3 . Pressure at collapse depth 
Young's modulus 

5 . Hull diameter 

6 . Conpartment length 

J. Density of the material 

8. Allovable out-of -roundness of the hull 



Page Two 



1 . Prame Area 

2. Number of frames 

3-8. Critical pressure predicted by formula 92a- 
modified^ Weak, Lunchick, and three modes 
by Bryant 

9-11. Three values of yield strength predicted by the 
formula in section l-c-6 of this appendix. 

12. Frame spacing 

13 . Shell thickness 

14. Web thickness 



15 • Volume of material in the shell and frames 
16 . Weight of material in the shell and frames 

17 - 38 . Preliminary calculation values beginning with "THETA" 

and ending with "C" as listed on the fourth page of the 
computer program. 

39 - 50 . Values used in the solution of the two differential 
equations when determining the required changes. 

51 . "Least weight" of those compartments considered. 



I 

I 



1 

I 



& 







d. Since thickness is first decreased to 95^ of the 
predicted value and then increased to of that 

predicted hy the hoop stress formula, the values on 
page two will he printed at least three times- If, 
as may happen in some cases, the solution with 
thickness increased produces the "least weight" 
solution, thickness will continue to he increased hy 
increments of k^o until a "least weight" solution is 
reached and passed, (Each of these incremented 
solutions will also require one page of printed 
output.) This series of events would occur if, hy 
increasing thickness, frame size and the number of 
frames could he decreased enough to achieve an 
overall reduction in weight. After a certain number 
of thickness increases, however, weight will begin 
to increase at which time the program will terminate. 
Experience indicates, however, that the "least 
wei^t" solution will occur when "T" is reduced to 
95^ of the thickness predicted hy the hoop stress 
formula, and if the hoop stress prediction is ignored, 
"least weight" will occxir at about 70 ^ of the 
predicted thickness. 
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YIcLC STRENGTH 


=: 


5 5 5 E . j E 


P.S. I . 


CCLL/lPSE DEPTH 


= 


153" . 


FEET 


PRESSURE AT CCLLAPSE DEPTH 


= 


68C . CO 


P.S. I . 


YCUNGS MCnULUS 


37 


0 OG C 0’'i . G'D 


P.S . I . 


f^EAN DIAYETER CF HULL 


= 


21A .70 


INCHES 


BULKHEAD SPACING 


= 


4 1 6 . : o' 


INCHES 


DENSITY CF NATERIAL 


= 


.2 340 


LB. PER 


ALLOWABLE CU T-CF-RCUNDN E S S 


= 


. 50 


INCHES 



- 6o - 



INCH 



FRAi^E AREA = 8.26 SG. IN. 

NC. FRAMES MET'aEFM BULKHEADS = 19.56 FRAMES 



CCLLAPSE7PRESS . FRUM 92AMCD. = 
CCLLAPSE DRESS. FRGM WENK 
COLLAPSE PRESS FROM LUNCHICK = 
CCLL PRESS FR BRYANT, 2 LC6ES= 
CCLL PRESS FR BRYANT, 3 LOBES= 
CCLL PRESS FR BRYAiNT, A LCBES = 
MAX STR FR 6USHIPS, 2 LCBES = 

MAX STR FR BUShIPS, 3 LCBES = 

MAX STR FR BUSHIPS, A LCBES = 

FRAME SPACING 
SHELL THICKNESS 
KEB THICKNESS 

VOLUME CF CNE COMPT MATERIAL = 
WEIGHT CF COMPARTMENT 
2 . A9059 
5.99276 
6.07562 
.81633 
1 .CA127 
1.2A530 
1.593j5 
1.68C91 
-.5:168 
-. 3A727 
. 19287 
.65292 
.1A018 
.C9028 
2.53025 
1 .COG 00 
. 85262 
. 89355 
. 6A815 
18. 18325 
131 . 12A29 
9.08A98 
■ A.7059A 

663, 1A37A 
55986. 18652 
-6. A2158 
3A7. 75732 
1.38360 
7. 5367A 
-219.65967 
-87.C5908 
1 lOA. 7082A 
3239.C3076 
2.93203 
6. 10A57 

COMPARTMENT CF LEAST WEIGHT 



l-'^j . '2 


P . S . I . 


853 . lA 


P.S.I . 


791 .7A 


P.S . I . 


A 7 2 3 . 3 1 


P.S.I. 


1558. 19 


P.S.I. 


2207 .23 


P.S.I. 


A 1390.5'* 


P.S.I. 


A6662.06 


P.S.I. 


5A0A2.2A 


P.S.I. 


22 .A1 


INCHES 


1.25 


INCHES 


.66 


INCHES 


228702.30 


CUBIC INCHES 


6A951 .A5 


POUNDS 



= 63257.29 POUNDS 
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